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Abstract. We extend earlier work on nonlinear tube wave propagation in permeable formations to study, 
analytically and numerically, the generation and propagation of a difference frequency, Auj = uji — uj2, due 
to an initial pulse consisting of carrier frequencies uji and uj2- Tube waves in permeable formations have 
very significant linear dispersion/attenuation, which is specifically addressed here. We find that the difference 
frequency is predicted to be rather easily measurable with existing techniques and could yield useful information 
about formation nonlinear properties. 



Introduction 

A tube wave is an acoustic normal mode in which the energy is confined to the vicinity of a fluid-filled 
cylinder within an elastic solid. From a practical point of view it is generally the dominant signal which appears 
in a typical borehole-logging measurement and thus it is important in a variety of contexts in the search for 
hydrocarbon sources. 

One of these contexts lies in the fact that the tube wave may couple to fluid flow within the rock formation 
if the latter is permeable. The linearized tube wave propagation in this regime has been extensively studied 



both theoretically and experimentally [Winkler, et ai, 1989] (see also [Pampuri, et ai, 1998] and references 



therein). In the present article, we use a model of the tube wave due to [ [Liu and Johnson, 1997| ]. According to 
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the model, the fluid in the borehole is separated from the porous formation by an elastic membrane (mudcake) of 
finite thickness. As a tube wave propagates, the membrane flexes in and out of the pores, thus forcing the fluid 
to flow through the formation. This leads to the coupling between the tube wave and the acoustic slow wave in 
the formation, which in turn leads to attenuation and dispersion of the tube waves. In formations of moderate to 
large permeability, this mechanism is the largest known source of attenuation/dispersion of the tube wave and is 
the reason why it is specifically considered in the present article. 

Quite apart from this effect it is also known that sedimentary rocks have very large coefficients of 



nonlinearity and so [Johnson, et ai, 1994] developed a theory for nonlinear tube waves neglecting the effects 



of the permeable formation. Later, [Johnson, 1999| combined this theory with the linearized theory of [Liu 



and Johnson, 1997] to describe a situation when the two effects are simultaneously present. As a numerical 



demonstration of the theory, [Johnson, 1999] considered the propagation of a narrow -banded (long duration) 
pulse consisting of a single carrier frequency. He showed that for realistic system parameters, the main signal 
(the fundamental) quickly decays, but before completely disappearing it generates a second harmonic and a 
low-frequency band (the "self-demodulated" pulse) both of which are due to the nonlinearity of the problem. The 
second harmonic decays even faster than the carrier, with the result that the self-demodulated pulse eventually 
dominates the entire signal at large enough distances. 

Because the second harmonic decays so fast, often it is advantageous to determine nonlinear characteristics 
by using pulses which consist of two different carrier frequencies, oji and uj2- In addition to the second 
harmonics (above) nonlinear effects lead to the generation of a signal centered around the difference frequency 
Auj = (jJi — UJ2- This component may be reasonably energetic while at the same time it is not attenuated as much 
as either second harmonic, or even as much as either carrier frequency. Thus, in this article we are motivated to 
consider the propagation of two narrow -banded pulses whose frequency separation Auj is, say, 10% of the central 
frequency. Moreover, because Aco is not that different from 2, it is often possible to measure its amplitude 
with the same acoustic transducers as for the fundamental. 

The organization of the article is as follows. First, we review the theory and derive analytical results for 
the nonlinear propagation of a pulse consisting of two different frequencies. We derive an approximation for 
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the propagation of the entire signal and we find an analytical form for the energy of the band with frequency 
Aw. Next, we report results of numerical calculations for a few different parameter sets and we show a good 
agreement between our analytical and our numerical results. In the last section, we make a brief summary of our 
work. 



Theory 



The dispersion and attenuation of the linear tube wave propagation has been studied in [Liu and Johnson, 



1997]. To simplify our discussion we use an approximate form of the dispersion relation from [ (ohnson, 1999] 



kz{uj) ^ UJ Sac + 0(w) 



(1) 



where is the wave vector along z-direction (the tube axis), lu is the angular frequency of the wave, Soo is the 
slowness at infinite frequency and 6 is given by 

Pf 



Soob[W^c + WpiL0)Y 

Here, p/ is the density of the borehole fluid, b is the borehole radius, Wmc is the mudcake membrane stiffness 



(2) 



defined in [Liu and Johnson, 1997], and Wp characterizes permeability effects. Wp depends on the borehole 
fluid viscosity, rj, formation permeability, k, and the diffusivity of the slow wave 

kK*, 



Cd 



(3) 



through the equation 



nHi'\ksib) 



Wpiu;) 



(4) 



Here, fcgi = \J ''■^ jCo is the wave vector of the slow-compressional wave and h'^^I are Hankel functions. 

Tube waves in nonlinear hyperelastic and impermeable formations have been studied in [Johnson, et ai, 



1994]. In [ lohnson, 1999| ] the effects of linear attenuation/dispersion and nonlinearity have been combined to 



obtain an approximate equation of motion for tube wave propagation in a realistic borehole. In the retarded time 
frame t = t — SoqZ this equation is 

dp{z,r) , dF{z,r) (]Si,dp^z,T) 



dz 



2pf dr 



0, 



(5) 
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where p is the pressure, is a dimensionless parameter defined in [Johnson, 1999 1. The function F{z,t) is most 
simply related to the acoustic pressure in the Fourier transform domain: 



Fiz,u;)^eiu;)piz,iu). 



After performing the Fourier transform of Eq. (|]), one obtains 



iq[u:) p{z, uj) + — (2, = 0, 

dz J 2pf 



(6) 



(7) 



where 

q{uj) = cje(w) (8) 

is the reduced wave number. In Eq. (0), p'^{z, uj) is the Fourier transform of the square of p{z, r) (not the square 
of the Fourier transform). 



In [ Johnson, 1999 |, Eq. (|7|) has been used to study the generation of the second harmonic as well as that of a 
low-frequency self-demodulated signal in a situation in which initially the pressure is a narrow band/long duration 
pulse centered on the frequency, ui. In the present article we consider a situation when two narrow-banded 
pulses are initially present. As discussed in the Introduction, for practical relevance we will take the difference in 
the pulse frequencies Auj to be 10% of the central frequency u;'. The initial signal is given by 



p{z = 0, t) = El (r) sin[a;iT + 0i] + £'2 (r) sin[ijj2T + (j)2] 



(9) 



Here, Ei{T) are the envelope functions, wi = oj' + Acj/2 and 0^2 = — Auj/2. For the numerical demonstration 



of the next section, we take the same envelope functions as in [Johnson, 1999|: 



Ei{t) = E2{t) = -Pocxp [-{r/Twr] ■ 



(10) 



By setting uji — uj2 and (pi ~ (j)2 in Eq. (M) one recovers the pulse considered in [Johnson, 1999]. In all three 



parameter sets used for the numerical calculations in the present article (see Table I), we take Auj/ui' = 0.1 and 
uj' X Tw — 1257r so that initially the different signals we consider all "look" the same. 

There are two characteristic distances relevant to the problem: the decay length of the linearized theory, 

Zatt = l/j{oj'), where 

7(w) = [e(w)j , (11) 
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and the distance over which a shock front would develop, in the absence of attenuation. The latter is found in 



[ Hamilton and Blackstock, 1998 1 to be 



^shock 



Pf 



The Gol'dberg number 



r — Ziitt/Zshock 



(12) 



(13) 



measures the importance of nonlinear effects relative to the linear. In the three parameter sets (Table 1), the 
amplitude Pq is chosen so that F = 0.21, i.e. the nonlinear effects are significant, but overall pulse propagation is 
dominated by linear dispersion/attenuation. 

For r -C 1, the nonlinear effects can be ignored to a first approximation. Because, by assumption, the signal 
consists of two narrow-band pulses one has the usual result of linear acoustics: 



p\z,t) = Ei{t ~ ASgiz)e-^'' sm[iJi{T - ASpiz) + <j)i] 

+E2{t - A5'g2^)e~'^=^'= sin [uo2{t - ASp2z) + 02] , 

where ASp (ASg) is the additional phase (group) slowness relative to 5*00: 

ASp{uj) = K \e){uj)] , 

d 



(14) 



ASgico) 



3? 



(15) 



Throughout this paper, the subscripts 1 and 2 correspond to frequencies uji and uj2, respectively. The physical 
meaning of Eq. ( p^ ) is transparent: The pulse envelope propagates with the group velocity and attenuates 
while each peak and trough travels with the phase velocity. If the envelope is broad enough, the relevant 
attenuation/dispersion quantities should be evaluated at the central frequency of the pulse. 

If the Gol'dberg number T is small compared to unity but is not completely negligible, Eq. ( [l4| ) still gives 
a reasonable approximation for the propagation of the original pulses but, because of the quadratic nonlinearity 
in Eq. (Q), new frequency components will be generated which are centered around uj = 2uji, 2uj2, as well as 
those centered around uji — U2 and uj\ + 102- Apart from a direct numerical solution of Eq.(^ one can develop a 
perturbation theory thereof: 
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p{z,t) =p'{z,t) +po(z,t) +Pa{z,t) +p2i^-,iz,T) +p2i^2{z,T) + p^-^+^2{z,t) + 



(16) 



Here, pq refers to the self-demodulated signal, pA refers to the signal whose bandwidth is centered around 
A — LOi — LU2, etc. By substitution of Eqs. ( [T4| ) and (|l6|) into Eq.(^ one can derive an approximate equation for 
the evolution of these nonUnearly generated components, even in the presence of significant linear dispersion and 
attenuation. 

In this manner the generation of the second harmonic and of the self-demodulated signal has been studied in 



[ [Johnson, 19990 . In ths present article we focus on the component centered around the frequency A = wi — 
(The analytical expressions we derive for this mode, can easily be generalized to the uj = uji + uj2 case by 
redefining lj2 —^2-) We obtain (At/) = — (/)2) 



Pa{z,t) = i?l(T - ASg^z)E2{T - ASg2z) ^'f^^^ X ^ 



g(Acj) -qi+q2 



(17) 



One is reminded that this equation is valid only for narrow -banded pulses (Auj x Tw ^ 1) and small nonlinear 
effects (r < 1). 

Next, we consider the energy of the carrier {uj « uj') band, £'{z) as well as that of the band centered on the 
difference frequency uj « Aco, £a{z)- We define them by 



£'{z) 

£Aiz) 



ZJ 12 

'72 
3A[^/2 

Aw/2 



\p(z^ui)\ dw, 



|p(z,w)l"'du;. 



(18) 
(19) 



Using Parseval's theorem and Eqs. ( 14 ) and (jlTp, one gets 

/ + 00 
\p'{z,u)fAu 
-oo 

/ + 00 /' + CXO 

Ef{T)dT + e-^^'' Ei{T)dT 
-oo J — OO 



(20) 



and 



£Aiz)/T: 



+ 00 



\pa{z,u)\ duj 



(AcvpSl 
V 2Pf 



^i(gi-gj)z _ pjq(Aw)z 



g(Aw) -qi+q*2 



Ef{T - ASgiz)Ei{T ^ ASg2z)dT. (21) 
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In the special case that the two envelope functions are identical, Ei{t) = E2{t), one can derive a simple 
analytic result for the nonlinearly generated signals centered around Aoj as well as the self-demodulated signal 
centered around cj = 0. For this, we rewrite the input signal in the form 

p{z = 0,T)^E'{T)s\TL[uj'T + (t)X (22) 

where = (0i + (?!)2)/2 and 

E'{t) = 2Ei{t) cos [(Acj/2)t + A0/2] . (23) 

E'{t) is now viewed as a narrow -banded envelope function for the w' mode (assuming that Aw ^ w')- As 
before, the total signal can be approximated by 

p(z,t) «p[5io^](z,r) +p'(2;,t) (24) 

where p[siow\ [z, t) = pq{z, t) + {z, r) now includes both the self-demodulated component as well as the 
components centered around Auj. Within the context of the foregoing approximations it is given by 

(3Sl /■+°° -tujE'^ioj) 



(25) 



Plsio^] {z, r) ^^^g, _ _ ^^^^^ 

The primed quantities AS'g and 7' ai-e evaluated at the frequency uj'. 

Numerical Results 

For the numerical calculations, we consider three parameter sets, {A, B, C}, which are listed in Table I. 
These parameters are identical to those considered in [ lohnson, 1999| , where the relevance of these parameters 



to realistic borehole properties is discussed. There is one change in that in this article we take Tw to be five times 
larger This is done in order to make the band widths of the pulses so narrow that the Aw mode clearly separates 
from the rest of the low frequency signal. 

In order to solve the equation of motion ([7|) we use the same Lax-Wendroff algorithm as in | ]johnson, 1999| ] 



and the initial pulse with envelope functions (|o|) and phases (f>i — <j)2 = 0. Figs, (jl]), (^, and (|]) show several 
snapshots of the pulse profile in both time and frequency domains for three parameter sets (Table I) A, B and C, 
respectively. Cases B and C have similar initial pulses, but differ somewhat in their dispersion relations, as can be 
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seen from the corresponding parameter sets in Table I. Sample A has a center frequency two orders of magnitude 
lower than samples B and C, and as a consequence, the relevant modes have much longer attenuation length Z^tt- 

It is pedagogically useful to first examine Fig. |]in detail. Considering the left column of plots, one can 
trace an intuitively clear sequence of events: Initially, the signal is given by two sharp pulses separated by 10% 
frequency difference Aui (first row of plots). As this signal propagates through the borehole and gradually 
attenuates, modes centered at w « Aut are being generated (second row). They soon start to dominate the general 
shape of the pulse because they attenuate less than does the carrier (third row). Concomitantly, w modes start 
to appear and give a significant "background" for the Acu signal (fourth row). In the end, all higher-frequency 
components decay and the completely self-demodulated signal w contains most of the energy (fifth row). 
One has to be careful looking at the right column of the plots: The black profile of the first plot outlines an 
envelope Eq. ( p3| ) of the lOkHz signal, detailed structure of which can be seen only if we expand the time scale r 
(second and third plots). In order to examine the Auj w lOOHz components and the lower-frequency signal, we 
again show the entire signal in the last two plots. Although Eq. ( |24| ) is also plotted, it is indistinguishable from 
the results of the full numerical calculation. For large enough distance p{z, t) evolves to P[siow] i^i Eq. 
as is indicated in the plot. 

Qualitatively Fig. |l] and Fig. ^ look similar to Fig. |^. But in the case of sample A, because of the smaller 
attenuation, the Auj component does not demodulate completely even after propagating a distance of ^ 40km. 

Next, for each of the three cases. A, B, and C, we calculate the energy of the ui — oj' band directly using 
Eq. (|l8]) and also using Eq. (^0|). Similarly, we calculate the energy of the uj = Aoj band using Eq. ( [l9| ) as well 
as Eq. (|T]). The results are plotted in Figs. ^ || and |[ One can compare Fig. ^ with Fig. 6 of [ lohnson, 1999 | 
which shows the energy of the second harmonic and the u; = band for an almost monochromatic initial pulse 
and same system parameters. As can be intuitively expected, the Alu band is intermediate in both the maximal 
energy it gains and the distance it propagates, as compared to the lower-frequency band and the higher harmonics. 
The Auj component does not reach energies as high as the second harmonics do, but it propagates much further, 
still carrying a significant fraction of energy; it is ^ 55dB down from the initial carrier energy, £'{z = 0), in our 
examples. On the other hand, while the lower-frequency signal persists longer than the Auj modes, it is never as 
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energetic and it would be more difficult to measure even in principle. 

The intersections of the solid and dashed curves in Figs. ^ ^ and ^ indicate the crossover region, where the 
Lu = Alu band starts dominating over the carrier band uj = lu', as can be seen in Figs. |l], ^ and|3], respectively. 

Conclusions 



We have used the theory of [ Johnson, 1999 | of tube wave propagation in permeable formations to describe 
nonlinear interaction of two narrow-banded pulses. The theory incorporates both nonlinear effects and a realistic 
model for dispersion/attenuation of tube waves. We have extended this previous work on the propagation of a 
single narrow-banded pulse to describe the generation and propagation of the lu — Aou band when two different 
carrier frequencies are present in the initial pulse. 

We have derived analytical results for the self-demodulated component, the Auj band, and the total signal in 
the regime of weak nonlinearity and they are in excellent agreement with an accurate numerical calculation using 
three different parameter sets. Specifically, we have studied the spectral content of the signal and demonstrated 
that the Alu band can have a potential application because of its long attenuation length and its relatively high 
energy content. Also, if Alu/uj' ^ 0.1, as we consider here, there is a practical bonus as the same transducers 
used for the generation of the carrier signal can, presumably, be used for the detection of the Aoj band. 
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Table 1. Values of input parameters for the calculation of tube wave 

characteristics. 



Sample 



B 



Input 



K* 



Pf 

Snc 



parameters /3 



u;'/27T 

Tw 

Po 



0.30 0.30 0.30 

0.01 0.01 0.01 poise 

2.25 2.25 2.25 GPa 

0.1 0.1 0.1 m 

1000. 1000. 1000. kg/m^ 

667. 667. 667. /isec/m 

50.5 50.5 50.5 

250. 250. 100. GPa/m 

2. 0.2 0.2 /Ltm^ 

0.1 10. 10. kHz 

625. 6.25 6.25 msec 



81. 



80. 92. kPa 



C 15. 1.5 1.5 m2/sec 

Zatt 278. 2.8 2.4 m 

Calculated Zghock 1323. 13.4 11.6 m 

quantities T = ^att/-^shocic 0.21 0.21 0.21 

ASp{uj') 55.9 6.8 6.9 /zsec/m 

ASg{uj') 52.7 3.5 3.5 /isec/m 
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Figure 1. Calculated pulse evolution using parameter set A, from Table I. The right column is the signal in the 

time domain, the left is the frequency spectrum. Also shown as a solid line on the right is the analytical expression 

for the low-frequency pulse, Eq. (^sf), toward which the pulse evolves. Similarly, the analytical expression for the 

total signal, Eq. (Q), is shown as a dotted line: It essentially overlies the numerically calculated signal. Notice 

the various changes of scale. 
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Figure 2. Calculated pulse evolution using parameter set B, from Table I. Same conventions as Fig. |l|. 
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Figure 4. Numerically calculated values of the energy £'{z) (solid line) and £a{z) (dashed line) defined by 
Eqs. and (|l|) for the waveforms of Fig. [l|, sample A. Each curve has been normalized by the value of 
£'{z ~ 0). The data are shown in both log-log and semilog plots, in order to emphasize the short and long 
distance behavior, respectively. The dotted lines show the analytical expressions, Eqs. ( ^0[ ) and (|2l|); the fact that 
they are almost indistinguishable from the numerically calculated curves is precisely the point. 
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Figure 5. Numerically calculated values of the energy £'{z) and £a{z) for the waveforms of Fig. ^, sample B. 
Same conventions as Fig. ^ 
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Figure 6. Numerically calculated values of the energy £'{z) and £a{z) for the waveforms of Fig. ^, sample C. 
Same conventions as Fig. ^ 



